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SUMMARY 


An account is given herein of an investigation conducted to deter- 
mine the feasibility of artificially thickening a turbulent bo\mdary 
layer on a flat plate. A description is given of several methods used 
to thicken artificially the boundary layer. It is shown that it is pos- 
sible to do substantial thickening and obtain a fully developed turbu- 
lent boxmdary layer, which is free from any distortions introduced by the 
thickening process, ajid, as such, is a suitable medium for fimdamental 
research. Measurements of mean velocity, spectrum of u-fluctuation, and 
intensity of u-fluctuation which served as criterions for determining 
the state of the layer are presented. Some features of the fully devel- 
oped turbulent boundary layer with zero pressure gradient are discussed. 
The mean flow in the boundary layer is compared with the logarithmic 
laws of mean- velocity distribution derived for pipes and channels, and 
the existing theories dealing with the shape of the spectrum in isotropic 
turbulence are applied to the spectrum, measured in the outer portion of 
the boundary layer. It is also shown that the turbulence in the outer 
portion of the boundary layer has an intermittent character similar to 
that found. for wakes and jets. 


INTRODUCTION 


The lack of a satisfactory theory for turbulent shear flow and the 
development of the hot-wire anemometer as a powerful tool for turbulence 
research have directed attention to the quantitative measurement of the 
characteristics of turbulence in shear flow as a means of obtaining 
information on which to base a sound theory. The recent concept of 
local isotropy proposed by Kolmogoroff (references 1 and 2 ); that is, 
that, whatever the nature of the large-scale motions, the small-scale 
tiurbulence tends to be isotropic, has intensified the experimental 
program with the hope that the existing theories of isotropic turbulence 
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may be useful. While a considerable amount of work has been done in 
various flow fields (references 3 , 4, 5, and 6 ), attention is here 
, centered on the boundary layer, not only because of its practical impor- 
,tance, but to provide information which can be compared with that obtained 
in other types of shear flow and thus show differences and likenesses 
and in particular show whether certain phenomena are universal. It is 
felt that the results of the measurements can best be compared with the 
theory of isotropic turbulence if such complicating factors as pressure 
gradient, curvature, surface roughness, and compressibility are elimi- 
nated. Furthermore, the turbulence should be fully developed and the 
Reynolds number should be sufficiently high to bring the domain of local 
isotropy within the range of measurement. In addition the boundary layer 
should be as thick as possible in order to facilitate the use of multi- 
wire hot-wire probes and minimize errors due to wire length. 

These requirements mean that the surface should be smooth and flat 
and that high Reynolds number should be obtained by thickness rather 
than by high speed. Such requirements are difficult to meet, except in 
wind tunnels especially constructed for the purpose. The use of one 
wall of a conventional wind tunnel does not appear to be a satisfactory 
solution because of the likelihood of extraneous effects. The problem 
of obtaining a thick boundary layer was encountered at the National 
Bureau of Standards when a turbiaent-boundary- layer investigation using 
a flat plate was, undertaken. The test section of the tunnel was 19 feet 
long, but it was not feasible to use a plate longer than 12 feet. It 
was felt that the need for a thicker boundary layer than one which would 
nomnally develop on 12 feet of surface was sufficiently urgent here, and 
possibly in other laboratories with even less length available, to warrant 
an investigation of the possibilities of artificially thickening the 
layer. This involved not only methods ,of thickening, but also sett in g 
up criterions by which to establish the identity of a fully developed 
layer over a smooth surface. 

In determining the feasibility and validity of ar\;ificial thickening 
a number of measurements had to be made, including mean velocities, inten- 
sity of u-fluctuations, and the spectrum of u-fluctuations. As a con- 
sequence it is possible to show how such quantities are affected by 
thickening devices and what values they assume when the layer is no 
longer influenced by its beginnings. It is obviously very necessary 
to know whether there exists a "flat-plate boundary layer" of laniversal 
character depending only on Reynolds number. This is shown to be the 
case insofar as the observed quantities provide a proper sampling. 

The investigation was conducted \mder the sponsorship and with 
the financial assistance of the National Advisory Committee for 
Aeronautics. 
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SYMBOLS 


X distance along surface from leading edge of flat plate 

y distance normal to sxirface, measured from surface of 

flat plate 

U mean velocity in boundary layer 

mean velopity in free stream. Just outside boundary layer 

u, V, w X, y, and 2 components of instantaneous tiirbulent- velocity 

fluctuations 


2 2 2 
’ U, V, w 

u' 

p 

V 

P 

Po 

‘^l 

Qo 

^s 


T 


mean-sq-uare values of u, v, and w 
root-mean- square value of u 
density of air 
kinematic viscosity of air 
static pressure 

static pressure at x = 10.5 feet 
free- stream dynamic pressure 

free-stream dynamic pressiare at x = 10,5 feet 

dynamic pressure in boundary layer as measured with a 
surface tube 

shearing stress 

shearing stress at surface 
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boundary- layer thickness 

boundary- layer displacement thickness 




nw 


■boimdary- layer momentum thickness 


Jo 


JL 

Ul 


(‘■a 


d.y 


Reynolds numher "based on boimdary- layer thickness 


Reynolds. number based on boundary- layer momentum 
thickness 


ra 


longitudinal space correlation coefficient of u-fluctuation 

longitudinal scale of u-fluctuation 

microscale in isotropic turbulence 

total dissipation of turbulent energy per unit mass 
per second 

time 

frequency, cycles per second 
one-dimensional wave number 
three-dimensional wave nimiber 

normalized one-dimensional spectral function associated 
with frequency n 
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Ex(kx) one- dimensional spectral function associated with 

E(k) three-dimensional spectral f\mction associated with k 


APPARATUS AND PROCEDURE 
General Arrangement 


The investigation was conducted in the 4i - f oot wind tunnel at the 

National Bureau of Standards, using a flat plate 12 feet long, moimted 
vertically along the center line of the test section, the leading edge 
being 6 feet from the upstream end of the test section. The plate was 
l/4-inch al\jminum with a smooth surface and a sharp leading edge formed 
by a smooth symmetrical ■ taper beginning 4 inches from the leading edge. 
The plate was fastened along the horizontal edges to one side of 3- inch 
channel irons anchored to the upper and lower tunnel walls, and so 
divided the test section into nearly equal passages. It was considered 
desirable to conduct the investigation in a tunnel of low turbulence so 
that the turbulence of the free- stream would be negligible compared with 
that in the boundary layers. The turbulence level of the tunnel was 
0.02 percent at 30 feet per second and 0.04 percent at' 100 feet per 
second. An elevation view of the tunnel is shown in figure 1. More 
details as to the arrangement of damping screens used in obtaining the 
low turbulence are given in reference 7* 

In order to obtain a condition of zero pressure gradient, a flexible 
second wall of sheet aluminum was mounted on the tunnel wall facing the 
working side of the plate. ■ It was attached to the ends of screws 
threaded into the tunnel wall, and thus provided an adjustable section 
to make the passage between the plate and the wall sufficiently diverging 
to offset the pressure drop resulting from boundary- layer growth. 

Another important consideration was proper angle of attack to 
afford smooth flow entry at the leading edge on the working side of the 
plate. This was accomplished by displacing the upstream edge of the 
plate and by partial blocking of the passage on the working side. Since 
the obstructions causing the blocking had to be placed farther downstream 
than any usable portion of the plate, the passage was extended 4 feet 
by joining the trailing edge of the plate with a section of plywood. 

Under these conditions the press\ire distribution along the plate 
was found to be as shown in f igiore 2( a) . From 2 feet on downstream 
the pressure is seen to.be uniform on the average with variations of 
about 1. percent of Qq about the mean. These variations could not be 
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removed by the adjustable wall because they were associated with a 
slight waviness of the surface of the plate. 

Figure 2(b) shows the treinsition region at a speed of 108 feet 
per second for smooth entry conditions at the leading edge. The begin- 
ning of treinsition corresponds to a Reynolds nijmber (based on length of 
laminar flow) of 3> 175,000, and here transition results from the growth 
of more or less regular laminar-boimdary- layer oscillations. 


Instriimentation and Techniques 

Traversing mechanism .- In order to traverse in a horizontal plane 
normal to the surface the unit shown in fig\ire 3 was used. A micrometer 
screw with 0.5-millimeter pitch is driven through a set of beveled gears 
and a shaft by a selsyn motor remotely controlled. To provide for ease 
of alinement, the pressure probes and hot-wire holders are attached so 
that they can move fore and aft, and rotate in a vertical plane. By 
means of the bracket shown in the figtire the unit is bolted to a stream- 
lined strut so that the end opposite the gears is forced against the 
s\irface. Channel irons anchored to the floor and ceiling of the tunnel 
served as guide rails for the strut, and traversing longitudinally was 
done by moving the strut fore and aft by means of a chain and sprocket 
manually operated from outside the timnel. 

Measurement of pressure distribution, transition, and mean 
velocity .- Measurements of pressiire distribution and transition were 
made with the pitot- static head shown in figure 4. The impact tube 
makes contact with the stirface, and the static tube is opposite the 
point of contact l/4 inch away from the sirrface. The static and impact 
tube are made of 0.04-inch nickel tubing with ,0.003- inch wall thickness. 
The static tube is made ■ according to the conventional design for such 
a tube, and the impact tube is flattened at the end to form a rectan- 
gular opening 0.0l4 inch wide. The static tube was used for the meas- 
urement of pressure distribution and both impact and static tubes for 
the measurement of transition. The head was attached by a long arm 
(containing the pressure leads) to the strut and could be moved back 
and forth in the x-direction while the forklike guides maintained proper 
relation to the siirface. The determination of the transition point with 
the surface tube depends on the increase in impact pressure close to 
the surface when the layer becomes turbulent. The abruptness of, the 
increase is shown in figure 2(b) . The beginning and ending of the 
transition region coincide with the beginning and ending of the rise. 

This was checked with the hot-wire head shown in figure 4. This sled- 
like head was attached to the strut so that it could be moved fore and 
aft like the pitot-static head. The hot-wire was at a fixed distance 
of 0.03 inch from the surface, and visual observation of the fluctuations 
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on an oscilloscope for different positions along the surface confirmed 
the indications of the surface tube. Measurements of mean velocity 
were made with a pitot-static combination similar to the one described 
above but without the guides. It is shown with the traversing imit in 
figxare 3* No correction was made for the effect of turbulence on the 
readings . 

Measurement of intensity and spectrum of u-fluctuation .- A descrip- 
tion of the method and the hot-wire equipment used for measuring the 
intensity of u-fluctuation is given in detail in reference 8. Platinum 
wire 0.0001 inch in diameter and 0.5 millimeter long was used. This was 
wire drawn by the Wollaston process, and the method of making a hot-wire 
probe was, to etch away the silver coating and then soft- solder the wire 
to the tips of fine sewing needles. The time constant of the wire was 
approximately 0.25 millisecond and was determined by the square-wave 
method described by Kovasznay (reference 9). Compensation for the time 
lag was made by a resistance-capacitance network in the amplifier. The 
response of the uncompensated amplifier was flat from 10 to 5000 cycles 
per second and fell off to about 50 percent at 10,000 cycles per second. 
The mean- square output was read on a microammeter working out of a 
thermocouple circuit. 

In order to measure the spectrum of u-fluctuation the output cir- 
cuit of the hot-wire amplifier was bypassed and the signal fed into an 
Electrical Research Products, Inc., sound analyzer with a frequency 
range from 10 to 10,000 cycles per second and a choice of fixed band 
widths of nominally 5 or 50 cycles per second. The rectifier type of 
output circuit of the analyzer was replaced by an additional stage of 
anqplifi cation and a thermocouple and meter. The purpose of this was 
to read mean- square currents and also to obtain a highly danrped indi- 
cator. In reading the output, averages were taken over a period of 
1 minute. The mean- square voltage corresponding to the compensated 
hot-wire signal passed by the band was obtained by first reading the 
hot-wire signal on the output meter and then matching this signal with 
a known sinusoidal signal fed into the system from, a General Radio 
microvolter. The known signal had the frequency corresponding to the 
center of the band. The mean-square hot-wire voltage per unit frequency 
was. then obtained by dividing by the effective band width. The effective 
band width was defined as the rectangular band width having the same 
area as the experimentally determined band shape. These were 

Effective Nominal 

(cps) (cps) 


7.07 

57.9 


5 

50 
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It was found to be advantageous to use the 7*07 band from 10 to 
1000 cycles per second and the 57*9 band from 1000 to 10^000 cycles 
per second. 


Artificial Thickening of Boimdary Layer 

It is at once obvious that objects protruding outward from a 
surface create more drag than the surface that they cover, and hence 
offer a means of producing a layer of de-energized flow in a compara- 
tively short distance in the streamwise direction. The energy of the 
main flow is diverted into eddy motion in wakes, and, while the flow is 
txxrbulent, it does not have the same dynamic structiire as a fiilly devel- 
oped turbulent boimdaiy layer over a smooth siorface. The question is, 
after leaving dra^-producing objects, what length of smooth surface is 
needed for the layer to become the same as a smooth- surface generated 
layer? It would seem that the length would be less than that required 
to generate a layer of the same thickness, but this is by no means a 
foregone conclusion. 

Since the ultimate aim of the work was to produce a layer approx- 
imately 3 inches thick, this thickness having been chosen as sufficient 
for the convenient use of hot-wire probes, no effort was made to find 
the maximum thickness that could be produced in a given length. Emphasis 
was placed rather on simplicity of method and adequacy of tests. 

The turbiilent boundary layer along the plate was first obtained by 
the conventional method of free transition, that is, from the instability 
of laminar-boimdary- layer oscillations. Then, various devices such as 
rods, screens, and distributed roughness were used on the laminar layer 
to force transition and de-energize the flow. Attention was given to 
maintaining the condition of zero pressure gradient for the modified 
layers by adjusting the "false" wall. The first attempt at thickening 
the layer was to use rods, which were 0.04, 0.11, 0.l8, and 0.25 inch 
in diameter. These were placed singly in the laminar layer and in con- 
tact with the svirface at k feet from the leading edge.. As will be seen, 
these were not satisfactory because of the length required for their 
distorting effect to disappear. The second method consisted of placing 
wire screens of various mesh sizes at different angles to the plate and 
close to the leading edge, so that the screen and plate formed a narrow 
wedge opening into the wind. A 24-mesh screen of 0.0075- inch wire 
diameter, which was 1 foot wide, making an angle of 0.5° with the plate, 
and extending from floor. to ceiling, was a satisfactory thickening 
device, but it was not considered practical because of the possibility 
of accidental shifting. The third and final method consisted of 
cementing sandpaper to the surface beginning at the leading edge and 
extending for distances up to 2 feet downstream. Various amounts of 
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sand roughness were tried until the layer was of the desired thickness. 
The final arrangement of the roughness was No. l6 floor sanding paper, 
cemented to the plate from the leading edge to 2 feet downstream, and 
extending from floor to ceiling. 

Initially the one -dimensional spectrum of u- fluctuations was thought 
to be the best over-all test to be applied to the resisting layers. 
Accordingly, for a basis of comparison, considerable attention was given 
to the spectrum of the fully developed layer resulting from free tran- 
sition. It was found later that the mean- velocity profile was a still 
more sensitive test of the condition of the layer. Finally, the dis- 
tribution of u' across the layer was added to the list of criterions. 
These three, compared in different parts of the layers and compared with 
other available data, served as the sample phenomena to determine the 
state of the layer. 


RESULTS AND DISCUSSION 
Mean-Velocity Measiirements 


With free transition beginning at the 5-foot position, as ' shown in 
fi^wre 2, mean— velocity distributxons were determined at x = T*5> Q.5 
9.5, and 10'.5 feet at a free-stream speed of 108 feet per second.^ These 
are shown in figure 5. A lower value of Uj_ would have shifted the 

transition point too close to the trailing edge of the plate. Figure 6 
shows the same velocity profiles on a nondimens ional basis. In order 
to avoid the arbitrariness of 6, the nondimensional distance from the 
wall y/6 is used. It is seen from figure 6 that downstream from the 
extreme position of the transition point, that is, x = 6.75 feet, the 
profiles exhibit similarity . The two-dimensional nature of the flow 
was checked by making velocity measurements 10 inches below and above 
the center line at x = 10.5 feet and comparing with the velocity, dis- 
tribution at the center line. Agreement among the three sets of distri- 
butions confirmed the two-dimensionality of the flow. 

The velocity measiu-emehts for the 0.04- and 0.25- inch rods are 
shown in figures 7 and 8. At a free-stream speed = 108 feet per sec- 
ond, the flow, after passing over the rods,. became ttirbulent in the 
separated layer and reattached to the surface, in the turbulent state. 

At a lower free-stream speed = 30 feet per second for the 0.04-inch 

rod, visual observation of an oscilloscope trace obtained by traversing 
downstream from the rod at a distance of 0.03 inch from the siirface, with 
a hot-wire anemometer, showed the reappearance of laminar-boundary- layer 
oscillations, indicating that, the separated layer returned to the surface 
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as a laminar layer without the occurrence of transition. This effect 
was also noticed hy Liepmann and Fila (reference 10). Figures 9 and 10 
give some of the representative profiles for the 0.04- and 0.25- inch 
rod in nondimensional form. The profile for x = 5.17 feet in figure 9 
differs from those at 6.0 and 10.5 feet hy more than the experimental 
scatter. This would place the beginning of similarity for the 0.04-inch 
rod between x = 5-17 feet and x = 6.0 feet. In figure 10 there is no 
evidence of the beginning of similarity downstream of the 0.25- inch rod, 
and comparison of the profile at x = 10.5 feet with other distributions, 
at the same Reynolds nuniber, which do exhibit similarity, indicated 
that, even at x = 10.5 feet, 6.5 feet downstream from the rod, the 
velocity distribution is not yet similar. Figures ll(a) and 11(b) show 
the distributions for the 0.04-, 0.11-, 0.l8-, and 0.25-inch rods at 
X = 6.0 feet end x = 10.5 feet, respectively. It is seen that the 
deviation among the profiles is less at x = 10.5 feet than at 6.0 feet, 
illustrating the trend toward similarity and the dependence of the simi- 
larity position on the rod diameter. The trend in figure 10 is in the 
opposite direction of a Reynolds number effect^ but the Reynolds nvunber 
change involved is small enough to be negligible compared with the dif- 
ferences indicated. It is expected that immediately downstream from 
the rod there would be an influence on dP/dx and dP/dy, but there 
is no evidence that it persisted to any great length downstream. If the 
similarity position for the 0.04-inch rod is taken at x = 5*5 feet, 
that is, 1.5 feet downstream from the rod, the ratio of similarity dis- 
tance to rod diameter is 450. Assuming a linear effect of the rod 
diameter on the similarity distance would put the similarity position 
for the 0.11-, 0.l8-, and 0.25-inch rods at x = 8.2, 10.9, and 
13.5 feet, respectively. This is applied only for the condition of the 
rod at 4 feet and = IO8 feet per second. The wind speed and the 

position of the rod in the layer, that is, the size of the rod cortpared 
with the boundary- layer thickness, may be other variables which have an 
effect. The conclusion from this is that too much length is required 
for distortions to disappear for rods to be practical thickening devices. 

The velocity distributions at a free- stream speed = 43 feet 

per second, for a 24-mesh screen, making an angle of 0.5° with the 
plate, are shown in figure 12. Although no summary plot is given, it 
may be stated that the: prof lies from x = 3*5 feet to x = 10.5 feet 
exhibited similarity. The boundary-layer thickness at x = 10.5 feet 
is approximately 3 inches as compared with 1.4 inches for the condition 
of free transition. The effect of wider angle' or greater solidity was 
to move the similarity position farther downstream and to give a thicker 
layer. For example, the profiles produced by a 24-mesh screen, at an 
angle of, 2.5° with the plate, did not show similarity at x = 5»0 feet. 

At X = 9*5 feet the flow was similar, and . the ' layer was 5*2 inches 
thick. Since no measurements were made between x = 5.0 feet and 
X = 9.5 feet, the position at. which similarity, first occurred is .unknown. 
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At an angle of 5° the flow was no longer similar for the length of the 
plate. The use of screens seemed to he an effective method of thicken- 
ing but was not further investigated because of the difficulties entailed 
in moimting a screen, keeping it under tension, and maintaining its shape 
for long periods of time. 

Figures 13, l4, and 15 show the velocity distributions downstream 
from the distributed roughness at free- stream speeds of 35, 55, and 
108 feet per second. Since the general characteristics for the three 
speeds are similar, the nondimensional profiles are given in figure l6 
only for = 55 feet per second. It is seen that a distance of about 
3.0 feet was required for similarity and that at x = 10.5 feet the 
layer was approximately 3 inches thick, which for = 55 feet per sec- 
ond is the thickness of a layer having its beginning about l4 feet 
upstream. In other words, the virtual origin of the layer was dis- 
placed 3.5 feet upstream of the leading edge. There was some gain in 
effective distance at the higher speed. At = 108 feet per second, 

the virtual origin was about 5.5 feet upstream of the leading edge. 

Comparison of the velocity distributions in the various artificially 
thickened layers is made with that of free transition, in figure 17, at 
approximately the same Reynolds number Rg. The agreement is excellent. 

To test the \iniversal character of the velocity distributions, it is of 
interest to compare them with those of other investigations. In fig- 
ure 18, the distributions in the layer thickened by sand roughness are 
compared with the measurements of Hama (reference ll) and Schultz-Grunow 
(reference 12). The two curves show the effect of Reynolds number, and 
the data are siifficiently consistent to conclude that the velocity dis- 
tributions obtained for the artificially thickened layer are character- 
istic of a turbulent boundary layer with zero pressure gradient and a 
smooth surface. 


Comparison with Logarithmic Laws 

Much of the theoretical work that has been done on the turbulent 
boundary layer is based on the application of the logarithmic laws of 
velocity distribution derived for flow in pipes and channels (refer- 
ence 13) . These involve the adoption of the mixing-length concept, a 
form for the shearing stress distribution, and the omission of the 
influence of viscosity on the turbulence. 

The wall-proximity law and the velocity-defect law for the velocity 
distribution are compared with experimental results in figures 19 amd 20. 
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The value of in the present measurements was determined from the 

measured velocity profile hy means of the Von Karmlin momentum equation 
which for zero pressure gradient gives 

de Tq 
dx 2q^ 

The meah-velocity distributions for the artificially thickened turbu- 
lent boundary layer with zero pressure gradient, together with those 
of Schultz-Grunow ( reference 12) and Freeman ( reference l4) , are plot- 
ted in the form of U/Ux against yU^/v in figure 19. In contrast 
with the measurements for 'pipes and channels’, it is seen that there is 
a systematic increase in the deviation from the logarithmic law as the 
Reynolds nvunber decreases. In the region where the law is valid fig- 
ure 19 gives values of 4.2 for A and 5*75 for B. The following table 
gives the comparison with the measirrements in smooth pipes and channels 
for wall proximity. 


Investigator 

Type 

. A 

B 

Nikiaradse (reference 15). 

Pipe 

5.0 

5.75 

Nikuradse (reference l6) 

Channel 

5.0 

5-75 

Donch (reference 17) 

Channel 

5.0 

5.75 

Laufer (reference 5) 

Chanriel 

7.2 

5.75 

Klebanoff eind Diehl 

Boaandary- layer 

4.2 

5.75 


The commonly used values for the experiments of Nikuradse in 
smooth pipes for points near the wall is 5.5 for B and 5.8 for A, 
but Millikan (reference l8) has shown that values of 5*75 and 5*0 are 
in closer agreement with the measurements. The broken line in figure 19 
shows "that the power law 





is a good approximation for the higher Reynolds number Rg: = 152,000 

throughout most of. the layers but it should be remarked that at the 
lower Reynolds mambers the agreement was not so good. The distribu- 
tions at lower Reynolds naambers, on a log- log basis, had a sinuous 
appearance which seemed to decrease with increasing Reynolds manibers. 
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In fig\ire 20, both boundary- layer and channel results are plotted 
in the form suggested by the velocity- defect law. It is seen that 
results for both the boimdary layer and the channel fit the straight 
lines up to y/5 = 0.15. The value of B seems to be universal for 
boundary layers and channels within this region. While results for 
pipes are omitte:d, it is remarked that the same holds true for pipe flow. 
It will be noted that the observed results are independent of Reynolds 
number when plotted in this way, indicating separate universal relation- 
ships of the form 


1 U 
Ut 


= f(y/5) 


A point to be made with emphasis is that boundary- layer profiles are 
not like channel and pipe flow profiles, and it is seen that an arti- 
ficially thickened boimdary layer conforms to the pattern established 
in figures 19 and 20 after the similarity stage has been reacfied. 


Spectrum and Intensity of u-Fluctuation 

It has been shown in the previous section that, at a sufficient 
distance downstream, the mean- velocity distributions assume a universal 
character. However, the lack of knowledge concerning the relationship 
between the turbulence and the mean flow prohibits the use of mean- 
velocity distributions as. direct evidence that the turbulence is fully 
developed. Since it is impractical to measure all the characteristics 
of the turbulence, and in view of the interest in measurements of the 
spectrum in shear flow, the one -dimens tonal energy spectrum was chosen 
as the basis of comparison for the turbulence in the artificially 
thickened layers with that resulting from free transition. 

The mean-square value of the. u-fluctuations was measured between 
frequencies- n and n + dn at values of n from 10 to 10,000 cycles 
per second. Over most of the frequency range the accuracy of the spec- 
ti*um measvurements is within ±10 percent. The accuracy was somewhat less 
at the low-frequency end below 100 cycles per second because of large- 
amplitude fluctuations and at the high end above 6000 cycles per second 
because of noise and possible wire-length effect. No wire-length cor- 
rections were made, and the noise level at 10,000 cycles per second was 
approximately 50 percent of 'the total signal. A direct over-all meas- 
urement of the intensity u'/u^^ was made during each spectrum determi- 
nation to compare with that calculated from the spectrum, and to serve 
as a check on the accuracy. In general, the agreement was within 
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5 percent. The spectrum curves are given in the nondimens ional form 
UF(n) / 

— against Here F(n) is the normalized spectinjm function 

defined by 



dn = 1 


and Ljj is the integral scale defined hy 



( 1 ) 


where R is the longitudinal space correlation. The value of L was 

X 

obtained from the intercept of the F(n) curve at n = 0 on the 
assumption that a space correlation was related to a time correlation 
by X = Ut and using the Fourier transform written in the following 
form: 



( 2 ) 


From equations (l) and (2), when n = 0> 




(3) 


The values of Ljj calculated by equation (3) for the turbulent boundary 
layer resulting from. free transition and the boundary layers artificially 
thickened with a 0.04-inch rod and by roughness as described earlier are 
given in figures 21, 22, 23, and 24. It is noted that, except near the 
wall. Lx is constant across the layer and is proportional to the 
boundary-layer thickness, having a value of L^^/b of approximately 0.4. 
It should be pointed out that the exactness of equation (2) is open to 
question at high turbulence levels, since it requires that the fluctua- 
tions at a point be the result of a pure translation with velocity U 
of a fixed turbulence pattern. 
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For the condition of free transition, figure 21 shows a set of 
curves at x = 10.5 feet, with = 108 feet per second, at different 

distances from the surface. It is seen that there is a systematic 
variation across the boundary layer up to y/6 = O.57. Since L was 

used as the reference length, all curves are brought together at zero 
frequency. The trend in going toward the surface is for the higher 
frequencies (smaller eddies) to have proportionately more of the energy. 

It may be that near the sirrface some of the smaller eddies can 
derive energy from the steep velocity gradient which exists and that, 
farther away from the plate where the velocity gradient is comparatively 
small, local isotropy begins at a lower wave number. It is hoped that 
future measurements of the spectrum of shearing stress will throw some 
light on this point. 

Comparison is made in figure 22, at x = 10.5 feet and y/b = 0.03 
and y/6 = 0.57, with the spectrums obtained in the boundary layer, 
resulting from free transition and those obtained in the artificially 
thickened layers using a 0.04-inch rod and sand roughness; The free- 
stream speed was 108 feet per second for the rod and free transition 
and 55 feet per second for the roughness so that the boundary-layer 
Reynolds number is approximately the same. The spectrums are similar 
to one another and the same family of spectrums exists in the artifi- 
cially thickened layers as for that resulting from free transition. 

Spectrum measurements made at various distances downstream from 
the O.Ok-inch rod, at y/6 values of 0.57 and O.k and = 108 feet 

per second, are given in figure 23- It will be noted that the curves 
for X = 6.0, 7*5/ and 10. 5 feet are similar. Although a single curve 
is drawn through the measured points, there is evidence of a Reynolds 
number effect in the downstream direction on the high-frequency, end of 
the spectrum. The energy associated with the higher frequencies is so 
3mall that it does not distort the normalized curve at the low-frequency 
end. The curve for x = k.17 feet at y/6 = O.k which is only 2 inches 
from the rod is definitely dissimilar. This curve is practically flat 
up to 1500 cycles per second. The same characteristic flatness was also 
observed at y/6 =0.03, so no particular significance should be attached 
"to the value of y/S = 0.^, as compared with the rest of the curves with 
y/6 = 0.57. Measurement of the spectrum for the condition of free tran- 
sition, at X = 6.9 feet, which was as close as it was possible to 
approach the transition, because of the back and forth movement of tran- 
sition, showed no evidence of the characteristic flatness observed close 
to the rod. In fact, it gave evidence that ; similarity of the spectrums 
began at the downstream extreme of the transition region. In the case 
of free transition, the spectrum was also measured in the laminar layer 
at X = k.5 feet. A very sharp peak was found at 280 cycles per second. 
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the predominant frequency of the laminar-houndary- layer oscillations, 
emd little energy was found at other frequencies. There was no evidence 
of this peak in the turbulent bovindary layer at x = 6.9 feet. It was 
thought At first that the flatness of the. spectrum observed near the 
rod was associated with the distorted mean- velocity profile, and that 
here was evidence of a connection between t\irbulence distortion and 
mean- velocity distortion. However, figure 2h, which pertains to thick- 
ening by roughness, shows no effect in the spectrum at the 3-foot posi- 
tion where the mean- velocity profile was definitely distorted. It would 
seem that the flatness of the spectrum measured near the rod is an effect 
peculiar to the rod and may be associated with eddy frequencies from 
the rod superimposed on a more or less normal spectrum. 

It may be expected that, in the turbulent boxmdary layer, at a 
disteince far enough away from the wall, sin extensive region of local- 
isotropy exists. It is of interest to see how well the spectrums in 
the artificially thickened txirbulent boiindaiy layer wnere the Reynolds 
number of turbulence is high agree with existing theories predicting 
the shape of the spectrum in isotropic turbulence. 

Before making the coiiparison it may be well to discuss, briefly 
some important featvires of the theory. A fundamental concept in the 
theories introduced by Kolmogorof f ( references 1 and 2 ) and Heisenberg 
(reference 19) is that energy passes into the spectrum by way of the 
larger eddies and is then passed on to smaller and smaller eddies. 
Appreciable dissipation by viscosity starts only when the eddies get 
small enough. Beyond a certain wave number, smaller eddies get their 
energy only from larger eddies and are supplied sttfficient energy to be 
kept in statistical equilibrivim even in cases where the larger eddies are 
decaying. When the, Reynolds number of turbulence u'L^fv is high, there 
will exist in the lower wave-number part of the equilibrium range a sub- 
range in which inertial forces predominate and the transfer of energy 
is equal to the total dissipation e, and at the high wave-number end, 
another subrange where viscous effects predominate. Heisenberg assumed 
that the transfer of energy at wave number k was caused by a turbu- 
lence friction produced by eddies with wave niimbers greater than k and 
that the total dissipation involved a turbulence viscosity plus ordinary 
viscosity in an expression of the’ following form: 


e = 




E(k") 
v”3 ■ 


dk" 


,2k'^(k0 dk' 
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Here the second term within the brackets represents the turbulence vis- 
cosity, and X is an absolute consteuit. The solution of equation (4) 
gives the spectnam in the equilibrium range the form 


E(k) 


= (i-r* 


-5/: 


_k 

ke 




(5) 


where kg is a wave number in the intermediate rauge where inertial 
forces and viscosity both play a role, and is given by 



The equilibrium spectrum given by equation (5) deals with a three- 
dimensional energy spectrum. In order to compare with experiment it is 
necessary to transform this equation to the one -dimensional form using 
the following relationship also derived by Heisenberg: 



The one-dimensional form of equation (5) is foimd to be 

(^) = xio ^1^ 



(7) 
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where 



and 

kgjj = 0.6h^ kg 


In the inertial range, that is, k^ « kg^, equation (7) hecomes 


Ex 



9 

110 


Al^x 


-5/3 


( 8 ) 


and in the viscous range, that' is, k^ » ks> equation (7) "becomes 


where 



(9) 


In figure 24, equations (8) and (9) are compared with the spectrums 
measured at different distances downstream from the roughness, at 
y/& =0.57 and = 55 feet per second. At the high- turbulence 

Reynolds niAmber of the bo\mdary layer, the existence of a -5/3 region 

“7 

is quite evident. The agreement with the n ' law should be taken with 
some reservation because of the lesser accuracy of the measurements at 
the high-frequency end of the spectmm. It should be mentioned in pass- 
ing that the slope of the spectrum curves at y/6 = 0.03, such as shown 
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in figure 22, is approximately one, where the -5/3 region exists at 
y/6 = 0.57. 

It is seen from figure 25 that in the intermediate range of fre- 
quencies eqiiation (7) does not agree so well as the interpolation 
formula 


Ex(kx) 


9 

110 




-5/3 


1 + 


^\S/3 

ksx/ 


-2 


( 10 ) 


suggested by Heisenberg. For convenience, equations (7) and (lO) are 
represented in figure 25 as 


and 



where 


(7a) 


(10a) 


F(z) = Mili 
^x 

and 


z 


nLx 

U 
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and Aq and Zg are experimentally determined. Formula (lOa), as 

shown in figure 26 , is also in good agreement with some xmpuhlished 
measurements made in isotropic txirbulence 50 inches downstream from a 

1-inch grid in the NBS 4^ -foot wind tunnel. The fact that the analyt- 
ical relation (7) does not agree with the measured spectrums so well as 
the empirical expression (lO) would indicate that equation (4) may not 
yet he quite exact. Expressions (7a) and (lOa) are identical in the 
extreme regions of the equilibrium range, and it is of interest to cal- 
culate the value of x from the constants Aq and Zg. Since 

Aq and Zg are determined from a nondimensional normalized represen- 
tation they are connected to X bjr the following relations. 

The term E(k) of the three-dimensional spectrum is defined such 

that 


Jo 


E(k) dk = 


u^ 




From relation (6), it can be shown that 


Ex(kx) dkx = i ^ E(k) dk 


pop 

so that, for isotropic tiirbulence where u = v = w , 


0“ 


0 


Ex(kx) ^x=\ 


and 




noo 


4 


F(z) dz = 


Jo 


uo 


Ex(^x) ^ 
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since 


4x 


2ffn 2ffz 
U ■ Lx 


^ ^ F(z) 

4 2n 



In the inertial range, 


F(z) 


AqZs 


5/3 -5/3 

z 


and 



5/3 






= Ai = 


110 

9 


■^o^s 


5/3 


— , . 5/3 




Sit^Lx/ 


( 11 ) 


It will Le noted that e and must be determined in addition 

to the quantities already mentioned before X can be calculated by 
equation (ll). For isotropic tiirbulence, the expression for the dis- 
sipation (reference 20) is 



( 12 ) 


where X is the microscale, obtainable from a spectrum measurement by 



(13) 


Partial, justification for the use of such relationships where the 
turbulence does not have general isotropy is afforded by figure 27, 
where it is seen that about 90 percent of the integral in equation (13) 
is contributed by frequencies above 500 cycles per second. This is in 
the region where the -5/3- and -7-power laws hold and is therefore the 
region of local isotropy. It is concluded therefore as a' reasonable 
approximation that X determined by equation (13) is the appropriate 
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and only scale entering into the viscous dissipation. However, much of 
the value of is contributed below 500 cycles per second, and the 

assumption that u^ + v^ + w^ equals as made in equation (l2), 

is not a siifficiently good approximation. Since v2 and ^ were not 
measured, the results of reference (6) were used to obtain 


u^ + v2 + w^ = 0.62 (su^) 

at y/6 = 0.57- The factor 0.62 is applied to the right side of equa- 
tion (12) and to equation (ll) where u^ appears. 

The value of X calculated from equation (ll) for the tiorbulent 
boundary layer artificially thickened by sand roughness, at x = 10.5 feet 
and y/S = 0.57> with the values of Aq and- Zg given in figure 26 , 
is 0.23 • The value of X calculated for the case of isotropic turbu- 
lence shown in figure 26 is 0. 65 . By a similar procedure X was cal- 
culated for the spectrum measured by Laufer in fully developed channel 
flow (reference 5)^ and a value of 0.52 was obtained. At present there 
are insirfficient experimental data to determine the reason for the var- 
iations obtained. It is of interest to mention that the turbulence 
u'Lx 

Reynolds numbers — - — were 1790,’ 550, and I90 for the boimdary layer, 

channel,- and isotropic turbulence, .respectively-, and that corresponding 
increasing values of X were obtained. ; , 

The insensitivity of the spectrum to distortions in the mean- 
velocity distributions as illustrated by figure 24 made it desirable 
to check further on the fully developed nature of the turbulence in the 
artificially thickened layer by measuring the distribution of u'/U^* 

Intensity measurements were made .at various distances downstream from 
the roughness at ' = I 08 feet per second, and these are given in fig- 

ure 28 . At X = 3-0 feet, 1 foot downstream from the edge of the rough- 
ness, the intensity is flat up to y/& = 0.4 and is apparently the 
result of the diffusion of the high turbulence generated by the rough- 
ness. From x=5-0 to x=10.5 feet, the distributions of u' 
exhibit similarity, which is the same region for which the mean-velocity 
distributions were similar. The change in turbulence level accompanying 
a change in speed from I 08 to 55 feet per second is also shown. It 
should be noted that there is no change in level accompanying the vari- 
ations in boundary- layer Reynolds number due to "changes in boundary- 
layer thickness. 
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An interesting feature of the turbulence in the boundary layer was 
its intermittent nature in the outer portions of the layer. This was 
also observed by Townsend for a wake (reference .21) and by Corrsin for 
a turbulent jet (reference 22). It seems 'to be characteristic of shear 
flow with a free boiuidary. Figure 29 shows two oscillograms taken at 
X = 10.5 feet, at y/& = 0.03 y/5 = 0.75. At y/5 = 0.03 there 

is no trace of intermittency, while at y/6 =, 0.75 it is quite marked. 
Observation of the oscilloscope trace indicated there is a gradual 
increase in the intermittency as the distance from the surface increased. 
It is planned in the future to make some quantitative measvirements of 
the intermittency, and its effect on the measvired turbulence intensity. 


CONCLUSIONS 


From tests made to determine the feasibility of artificially 
thickening a turbulent boundary layer on a flat plate, the following 
conclusions can be drawn: 

1. It has been established that by artificial methods it is pos- 
sible to thicken substantially a fully developed turbulent boundary 
layer on a flat plate. The layer so formed becomes characteristic of 
that developed from the beginning on a smooth flat plate after passing 
over the required length of smooth surface. For a given boundary- layer 
thickness, artificial thickening makes possible a saving in length 
provided that there is enough length for the layer to become free from 
distortions. 

2. Almost any kind of obstruction placed on the surface would 
produce a thickened turbulent flow, but only from a certain class did 
downstream effects disappear in a reasonably short distance. This class 
appeared to be characterized by the ability to absorb energy from the 
mean flow without introducing large-scale disturbance. Two examples of 
this class are as follows: . 

(1) A fine-mesh screen placed so that the plate and screen 

form a narrow wedge opening into the wind 

(2) Sand roughness cemented to the surface 

3. The mean- velocity distribution appears to afford a simple means 
of identifying a fully developed turbulent boxandary layer. 

4. The logarithmic laws of velocity distribution are not in good 
agreement with the measured distributions in the t\irbulent boundary 
layer with zero pressure gradient and are applicable only in a small 
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region near the wall. The velocity profiles in the boundary layer are 
independent of Reynolds number on a velocity- defect basis, indicating 
a universal relationship of the form 


where U is the mean velocity in the boundary layer, is the mean 

velocity in the free stream just outside the boundary layer, is 

the square root of the ratio of the shearing stress at the surface to 
the density, y is the distance normal to the surface measured from 
the surface of a flat plate, and 6 is the boundary- layer thickness. 

5. The longitudinal scale of turbulence in a fully developed 

layer is proportional to the boundary- layer thickness 6, with a value 
of Ljj/S of approximately 0.4, and except near the wall is constant 

across the layer. 

6. There is a systematic variation in the one-dimensional energy 
spectrum across the boundary layer up to y/6 = 0.6, and in the outer 
region of the boundary layer, that is, y/6 >0.6, the existing theories 
for the spectrum in isotropic turbulence can be applied. 

7. There is an intermittent character to the flow in the turbulent 
boTondary layer similar to that found for wakes and jets. 


National Bureau of Standards 

Washington, D. C., September 15> 1950 
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Figure 3 .“ Traversing mechanism with total-head and static tubes attached 

Reference scale in inches. 
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Figoire U.- Sliding hot-vd.re head and surface pitot-static tube 
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Figure 5»“ Mean-rvelocity distributions for free transition as shown in 
figure 2, U]_ = 108 feet per second. 
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Figure 6.- Nondimensional mean-velocity distributions for free transition 
as shown in figure 2., = 108 feet per second. 
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Figiore 7 .- Mean-velocity distributions for O.Olt-inch cylindrical rod in 
contact with surface U feet from leading edge, U]_ = 108 feet per 
second. 
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Figure 8.- Mean-velocity distributions for 0.25-inch cylindrical rod in 
contact with surface U feet from leading edge. U]_ = 108 feet per 
second. 




Figxire 9.- Nondimensional mean-velocity distributions for O.OU-inch 
cylindrical rod in contact with surface 4 feet from leading edge 
Ui = 108 feet per second. 
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Figure 10,- Nondimensional mean-velocity distributions for 0.25-inch 
cylindrical rod in contact with surface k feet from leading edge. 
Ul = 108 feet per second. 
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Figure 12,- Mean-velocity distributions for 1-foot span of 2li-niesh 
screen making angle of 0.5° with plate. = U3 feet per second. 














Figure l^.- Mean-veloci.ty distributions- for No. 16 sandpaper covering 
first 2 feet of surface, U]_ = 108 feet per se.cond. 





Figure 1?.- Comparison of nondimensional mean-velocity distributions having 

approximately same values of R©. 
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Figure 19 . - . Coitparison of experimental results, with l/ 7 -power law' 

logarithmic law. 
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Figvire 20,- Channel and boundary-layer profiles compared on basis of 

velocity-defect law. 
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Figure 21,- Spectrum of u-fluctuations in turbulent boundary layer 
resulting from free transition, x = 10.5 feetj = 108 feet per 
second. 
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Figure 22.- Spectrum of u-fluctuations in turbulent boiandary layer at 
X = 10.5 feet. Ui =• 108 feet per second for free transition and 
O.OU-inch rod; Up = 55 feet par second for roughness. 
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Figiare23,- Spectrum of ; u-fluctuations in turbulent boundary layer for 
0,04-inch rod .in contact /with surface ’4 feet- from leading edge; 

U]_ = 108 feet per second. ■ ’ ■ ' ' ' ■ 
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Figure 2I4.- Test of -5/3 law and -7 law for spectrum of u-fluctuations 
in turbulent boundary layer thickened by roughness for first 2 feet 
of surface, = 55 feet per second. 
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Figvire 25.- Comparison of measured spectrum of u-fluctuations with 
Heisenberg's theoretical formula (broken line) and interpolation 
formula (solid line). Boundary layer thickened by No., 16 sandpaper 
for first 2 feet of surface, x = 10.5 feetj Ui = 55 feet per second. 
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Figure 26,- Spectrum of isotropic turbulence compared vjith Heisenberg's 
interpolation formula. Measurements made ^0 inches from 1-inch grid 
in -uniform stream, ' 




by formula 
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Figure 29.- Oscillograms of u-fluctuations in turbulent boundary layer. 

X = 10.5 feet; = 55 feet per second; 6 = 3.1 inches. Time interval 

between dots, l/60 second. Boundary layer thickened by No. 16 sandpaper 
over first 2 feet of surface. 



